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Coarse-grained spin density functional theory (SDFT) is a version of SDFT which works with number/spin 
densities specified to a hmited resolution — averages over cells of a regular spatial partition — and external 
potentials constant on the cells. This coarse-grained setting facilitates a rigorous investigation of the math- 
ematical foundations which goes well beyond what is currently possible in the conventional formualation. 
Problems of existence, uniqueness and regularity of representing potentials in the coarse-grained SDFT set- 
ting are here studied using techniques of (Robinsonian) nonstandard analysis. Every density which is nowhere 
spin-saturated is V-representable, and the set of representing potentials is the functional derivative, in an ap- 
propriate generalized sense, of the Lieb interal energy functional. Quasi-continuity and closure properties 
of the set-valued representing potentials map are also established. The extent of possible non-uniqueness is 
similar to that found in non-rigorous studies of the conventional theory, namely non-uniqueness can occur for 
states of collinear magnetization which are eigenstates of S^. 
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1. INTRODUCTION 

Modern electronic density functional theoryir— (DFT) 
is a very successful basis for the computation of many 
ground state properties used by chemists, physicists and 
materials scientists. However, it is not just practical com- 
putational algorithms, but also a distinctive way of look- 
ing at quantum many-body systems with its own set of 
fascinating basic questions. 

While the computational side of DFT is highly devel- 
oped, the mathematical foundations are relatively im- 
poverished. The fundamental questions revolve around 
the concept of y-representability. For a system of 
N electrons, with their mutual Coulomb repulsion, a 
single-particle density p{x) is said to be pure-state V- 
representable if some external single-particle potential 
v{x) has a ground state wavefunction with single-particle 
density p. It is mixed-state V-representable (simply "V- 
rcprcscntable" here) if there is a mixed ground state (den- 
sity matrix) with density p. Natural questions are those 
of existence, uniqueness and regularity of the represent- 
ing potentials as a function of density. The realization 
that some densities are not pure-state V-representable, 
but are mixed-state V-representable was responsible for 
the rise to prominence of the latter concep1j2r2i around 
1980. About the same time, some densities which are 
not even mixed-state V-representable were suggestedi^. 
Although V-representable densities are dense in an ap- 
propriate topology, so are non- V-representable densities 
and there is as yet no nontrivial positive characterization 
of V-representability. Concerning uniqueness, the origi- 
nal paper of Hohenberg and Kohni^ already contained a 
strong uniqueness result (modulo a trivial constant shift 
of the potential). Subsequently, when the theory was 
generalized to spin density functional theory (SDFT) in 



order to study magnetic phenomena, it was realized^^ 
that there was some breakdown of uniqueness in that 
context. There has recently been clarification^^"— of the 
extent of that nonuniqueness. Questions of regularity 
— how much will the representing potential of a density 
resemble those of nearby densities — has been almost en- 
tirely neglected. And not just by a lack of results; even 
the question seems largely unacknowledged. These foun- 
dational issues arc also relevant to the computational side 
of DFT. Computability is threatened not only by lack of 
existence, but also by lack of regularity. 

Chayes, Chayes and Ruskaii^ studied DFT for a lat- 
tice version of quantum mechanics and showed that, al- 
though serious mathematical difficulties or even patholo- 
gies may arise from the presence of arbitrarily short dis- 
tance scales, the infinite- volume limit is a tractable prob- 
lem. In particular, it was shown that every density is 
a ground-state density of an essentially unique poten- 
tial. On a lattice, short distance scales are completely 
eliminated; this is a fundamentally different quantum 
mechanics than the orthodox continuum version. An- 
other way to keep short-distance-scale degrees of free- 
dom from causing problems without altering the under- 
lying continuum quantum mechanics was introduced as 
coarse-grained DFT—. The idea is that we only allow 
ourselves to specify densities with some limited spatial 
resolution. On scales finer than the resolution, the den- 
sity is automatically relaxed by an energetic criterion. 
This approach seems much in the spirit of DFT. Con- 
ventional density functional theory asks for the lowest- 
intrinsic-energy (kinetic plus Coulomb) state consistent 
with a fully and exactly specified density. But it is en- 
tirely natural, and perhaps more computationally rele- 
vant, to consider incomplete specifications. Unique V- 
represcntability holds also in this formulation and fur- 
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thcr regularity results were recent l y^^i^^ demonstrated. 
Although the infinite-volume limit was handled, it was 
not found to be trivial, and for each property had to 
be approached anew. Since practical computations work 
with limited representational resources, usually in a form 
that amounts to limited spatial resolution, the coarse- 
grained formulation might be considered a more suitable 
grounding than the conventional fine-grained one. On 
the other hand, the coarse-graining scale can be taken 
as small as one wishes (10~^° m, say), so the approach 
is not inherently imprecise. Still, the infinite-resolution 
limit in which the coarse-graining scale goes to zero is of 
interest to understand the scale-dependences and to make 
contact with things formulated in a conventional infinite- 
resolution form such as an exact Coulomb potential. For 
that limit the coarse-grained formulation seems to have 
an advantage over the lattice, since a finer scale on the 
lattice requires a new lattice, not just some extra resolu- 
tion. Some progress in that vein has been mad o^^'^^ for 
ordinary (non-spin) DFT. This paper has no new contri- 
bution concerning that problem, but instead returns to 
the single-scale setting and the infinite- volume limit for 
spin-density functional theory. 

In previous work on the coarse-grained approach, it 
was observed that for a system confined to a finite box 
the situation really is simple. This paper aims to eas- 
ily obtain the infinite- volume limit by moving from finite 
boxes to the intermediate stage of a hyperfinite box. "Hy- 
perfinite" is meant here in the sense of nonstandard anal- 
ysis (NSA)^"— which is a rigorous way to use the notion 
of infinitesimals as well as their reciprocals. In that set- 
ting, a hyperfinite box can be asserted meaningfully to 
be larger than the infinite- volume limit, yet formal prop- 
erties of a system confined to such a box can be deduced 
immediately from those in a finite box. The added clar- 
ity brought by the infinitesimal methods allows not just 
ordinary DFT, but also spin density functional theory to 
be handled. Everything is done from scratch; there is 
no dependence on the earlier results. Thus, this paper 
also aims to promote the use of infinitesimal methods, 
which are not nearly as well-known as they deserve. It is 
hoped that nonstandard analysis will in the future allow 
the continuum limit also to be dealt with. 

The next section reviews some basic ideas of SDFT and 
Section|3]is a very brief introduction to nonstandard anal- 
ysis (NSA). Readers with appropriate background can 
skip these sections. Section |4] gets down to work, prov- 
ing some fundamental lower semicontinuity and near- 
standardness results using infinitesimal tools. The basic 
ideas and notations for the coarse-grained formulation 
are given in the short Section [S] 

The main results of the paper are found in Sees. IHl 
[7] and m The coarse-grained versions of the Lieb inter- 
nal energy functional F[p] and the set-valued function 
Potl (p) which gives the representing potentials are the 
main objects of interest. Section [5] proves that F is con- 
tinuous, that Potl (p) is non-empty if p is everywhere 
non-zero and nowhere spin-saturated, that Potl (p) is the 



functional derivative of in a sense appropriate to a gen- 
eral convex functional, that the map p i—^ {p ■ v : v € 
Potl(p)} restricted to nowhere spin-saturated densities 
is an upper semicontinuous set-valued function and 
that the graph of Potl is closed. The reason we have 
to work with Potl as a set-valued function is the well- 
known non-uniqueness in the SDFT context. Sections 
[7] and [S] deal with this problem. There are no nonstan- 
dard arguments in these sections, so they could probably 
be read on their own. The extent of non-uniqueness in 
the coarse-grained theory is shown to coincide with ear- 
lier non-rigorous conclusionai^ for the continuum theory. 
Namely, number/spin potentials may be non- unique only 
if the spin-density is saturated somewhere, or in case of 
coUinear magnetization in an eigenstate of Sz. The ap- 
propriate conditions are formulated in terms of densities 
rather than wavefunctions. Some concluding remarks are 
found in Section [5] Sections [7] and [S] have a very differ- 
ent fiavor from that of Sections 2] and IHl but they are all 
important parts of a well-rounded picture. 

The reader who is curious about the nonstandard ar- 
guments might consider reading quickly through Section 
[3] and then skipping to Sections [5] and [51 The reader who 
just wants to see the results may wish to start with the 
summary at the beginning of Section [SJ working back- 
wards as needed. 



2. SOME BASIC IDEAS OF SDFT 

In this section, we review some basic ideas of non- 
relativistic Density Functional Theory. It can safely be 
skipped by anyone with an acquaintance with that for- 
malism, after taking note of our notation for spin densi- 
ties. The discussion is kept at a heuristic level. 

We are concerned with a system of N identical parti- 
cles interacting with each other and subject to an exter- 
nal single-particle potential which functions as a control 
parameter. The wavefunction of a pure state for this sys- 
tem is a function ip{zi, Z2, . . . , z^) of N positions Xa and 
N spin components Sa with respect to some quantiza- 
tion axis. These are combined in the abbreviated nota- 
). In the usual concrete situations, the 
particles are electrons interacting via Coulomb repulsion. 
These are spin-1/2 fermions, so that the wavefunctions 
are required to be antisymmetric under interchange of 
and z/3 for a ^ p. The value of the particles' spin, and 
even whether they are fermions or bosons plays no cru- 
cial role. The discussion in this paper is tailored to the 
spin-1/2 fermion case, but appropriate modifications can 
be made for others. With the inner product 

= J ^{z)*4>{z)dz^---dzN, (1) 

where J dz denotes integration over position and summa- 
tion over spin, the antisymmetric wavefunctions comprise 
an Hilbert space which will be denoted i^. 
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As alluded to in the introduction, a general Density 
Functional Theory working only with pure states does 
not get very far. For reasons of convexityii it is advanta- 
geous to allow mixed states, also called density matrices; 
this is also a physically reasonable extension. Mathemat- 
ically, a mixed state is represented by a positive trace 
class operator in f), having an eigenfunction expansion 
of the form 

oo 

7 = ^c,|^',)(V',;|, (2) 

■i=i 

in Dirac notation, where tpi is an orthonormal set. 

The (single-particle) number/spin density correspond- 
ing to pure state 4' is 

Pap{x) ^ N J dz2- --dzN 

Ipixa, 2:2, ... , ZNTlpixP, Z2, . . . , 2Ar). (3) 

This is a 2 X 2 matrix- valued function of position. With 
(To the 2x2 unit matrix and ai, i = 1,2,3 the Pauli 
matrices, the number/spin density can be written as 

P = Pq<Jq + Picri + P20-2 + P30-3 = ^ ^Tr {a^p)cri. (4) 

It is therefore convenient to express the number/spin den- 
sity as the four- vector 

p := (po,p) = (po,pi,P2,P3), (5) 

and refer to it as the ^-density^ the first component being 
number density and last three proportional to the spin 
density. An alternative notation (n, 2m) = {po,p) will 
also be used, mostly in Sections [7] and [5] The factor of 2 
means that m actually integrates to the spin in units of h. 
The 4-density satisfies \p\ < po/2 everywhere. Notation 
like "ip I— >■ p" is customarily used to indicate that ip gives 
the 4-density p, but this is sometimes awkward, so we 
give the map a name. The 4-density corresponding to 
ip is denoted Dens-!/;. This statc-to-density map extends 
additively to mixed states; that is, 

Dens 7 = Ci Dens t/i,. 

i 

A physically normalized pure state satisfies \\ip\\ ~ 1, 
so that J podx = A^. Constantly making sure of the 
normalization is a distracting and unncccessary nuisance. 
Thus, we do not generally assume or insist that states 
and densities be physically normalized, unless otherwise 
noted. 

It is traditional to work with the ingredients of the 
Hamiltonian as Hilbert space operators, but quadratic 
forms are mathematically convenient and arguably 
more physically meaningful. Define the kinetic energy 
quadratic form by 

S'^ii^, VV'(^)* • dz, (6) 



the Coulomb interaction energy quadratic form by 

S^i^P,^) -.^Y. I \ M dz, (7) 

and the total intrinsic energy by 

£:(V',</>) :=f^(V,'^)+f'^(V','^). (8) 

For the moment, we ignore domain questions. The ab- 
breviation £{il),4) = £{'4') is used for diagonal elements, 
and similarly for £^ and f 

The Lieb internal energy functional^rii for SDFT is 
defined by 

:= inf{£(7) : 7 G 6,Dens7 = p}. (9) 

This is the minimum intrinsic energy consistent with den- 
sity p. The Lieb functional is simple and natural in 
retrospect, but took a long time to emerge. It solved 
the original "V-representability problem" , by extending 
the Hohenberg-Kohn internal energy functional to densi- 
ties regardless of V-representability (indeed, without any 
mention of potentials). To discuss the continuity of F, 
this most central object of DFT, requires a topology. One 
which suggests itself is the norm topology which is im- 
plicit in the very concept of density. In other words, we 
view the 4-densities as a subset of X : = L^(R^; C*) with 
the norm 

11/11 = Jmx)\ + \fix)\)dx. (10) 

Then, F is extended to all of X with value -l-oo off the 
range of the Dens map; this is just a convenience with no 
physical significance. There are highly oscillatory densi- 
ties in the range of Dens for which F = -l-oo, as well. 

Because the Lieb functional is defined in terms of 
mixed states rather than pure states, it is easy to see 
that F is convex: 

F[Xp + (1 - A)p'] < XF[p] + (1 - A)F[p'], (11) 

for < A < 1. This is an extremely important property 
and a nice way to rephrase it is in terms of the epigraph 
of F, 

epi^^:={(p',y):y>^^[p']}, 

which is the region on or above the graph of F in X x R. 
Thus, epiF is a convex set. It is also closed in the 
topologyii (see Cor. 14. 2^ : this property corresponds to 
lower semicontinuity of F. 

A representing potential v for p is a four-component 
function defined by the property 

F[p] + {v,p)<F[p'] + {v,p') yp'eX. (12) 

The notation {v, p) means 

(t!,p):= v pdx= {vopo + V ■ p) dx. (13) 
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Physically, p may be thought of as a magnetic field, 
though its divergence is unconstrained. 

It appears at first sight that some abstract results of 
convex analysis^ can be applied to the situation. For 
example, the Hahn-Banach theorem asserts the existence 
of a closed hyperplane separating (p, from the in- 

terior of epi F. Such a hyperplane would be the geo- 
metrical counterpart of a representing potential for p 
in X*, the dual space of X. Similarly, there is a the- 
orem (Prop. 1.5.3 o^) which says that when there is 
a unique such hyperplane at {p,F[p]), then it is the 
Gateaux derivative of F at p, if F is continuous at p. 
Unfortunately, these theorems get no traction whatever 
because F is not continuous. In fact, the effective do- 
main of F, domF = {p & X : F[p] < +00} has empty 
interior. This is easy to see^ii^. Any open ball around 
any point in X contains elements of X which are nega- 
tive somewhere, hence not even in dom.F. In addition, 
adding short-wavelength wiggles to a density can drive 
its kinetic energy arbitrarily high, or even to infinity, 
with arbitrarily small change in X-norm. Lieb^i took 
n L^{R^) for X instead of ^^(M^). does not 
contain the full range of Dens, but it does contain the 
densities of finite internal energy, so is acceptable. This 
choice of X, however, does not solve the problems just 
mentioned. Furthermore, there are reasonable potentials 
which are not even linear functionals. For instance, a har- 
monic potential \x\'^ takes value +00 on some densities. 
Eq. ([T^ makes perfect sense for such a case. 

The only general result on V-representability in the 
conventional framework is the fac t^^'^^ that the set of 
densities which are representable by a potential in X* is 
dense in domi^ in X-norm. This is a very weak result, 
and does not appear to be useful. 

Given a pair p, v satisfying Eq. ([T^ , v can be regarded 
as some sort of "derivative" of F at p. But what sort? 
In particular, we would like to know whether {—v, 5p) 
coincides with the directional derivatives 

F'[p; Sp] := hm - {F[p + edp] - F[p]) . (14) 

There is no general guarantee of this, as suggested above 
in reference to Gateaux differentiability. 

The questions just raised — which densities are V- 
representable (not to mention, just what qualifies as a 
potential)? to what extent do representing potentials co- 
incide with directional derivatives of F7 — are basic el- 
ements of an inquiry into the characteristics, potentially 
pathological, of the internal energy functional F. At least 
insofar as computational algorithms applied to practical 
problems are iterative, these questions are relevant, since 
the behavior of F in a full neighborhood of a sought den- 
sity is so. The rest of this paper aims at, and obtains, sat- 
isfactory answers to these questions in a coarse-grained 
framework, which is slightly different and less demand- 
ing than the conventional fine-grained one discussed in 
this section. Formally, everything is very similar, as the 
changes amount to a rcintcrpretation of density and po- 



tential. With a little nonstandard analysis, the answers 
for the coarse-grained theory are obtained easily. The 
next section therefore gives a whirlwind tour of the ba- 
sics of nonstandard analysis. 



3. SOME BASIC IDEAS OF NSA 

In the 1960's Abraham Robinson^^ revived the old 
idea of infinitesimals using methods from model theory 
(a branch of mathematical logic), and thereby creating 
the field of Nonstandard Analysis or infinitesimal analy- 
sis. Although they are a characteristic feature, however, 
infinitesimals are far from the only "ideal" objects Non- 
standard Analysis offers for dealing with mathematical 
problems. In this section, I try to prepare the reader with 
just enough of the jargon and basic ideas that the flavor 
of arguments in the rest of the paper may be appreciated. 
This flavor — once one gets used to it — is highly intu- 
itive. For further background, I highly recommend the 
notes of Lindstroro^i. For a highly compressed (twelve- 
page) summary, see Section 1 o^. There are also several 
good books2^"— , as well as some resources^^iM available 
on the web. 

In (model theoretic) nonstandard analysis, there is a 
nonstandard counterpart of every conventional mathe- 
matical object. A natural place to start is with the tra- 
ditional reals, M. Its counterpart is *M., the ordered fleld 
of hyperreals, which is an extension of R: in addition to 
the familiar real numbers, *R contains infinitesimals and 
their reciprocals, the illimited hyperreals. e S "^K is in- 
finitesimal if |e| < 6 for every (5 > in R, whereas x is 
illimited (called "infinite" by many) if |a:;| > n for every 
natural number n. The notation "n « -foo" (n « — cxd) 
means that n is positive (negative) illimited, and x 00 
means that x is limited. 

Two hyperreals x and y are infinitely close, x ~ y, if 
their difference is an inflnitesimal. Thus, "a; is infinitesi- 
mal" and "x w 0" are synonymous, x ^ y means that x 
does not exceed y by more than an infinitesimal. 

In fact, every object (formally viewed as a set) of con- 
ventional mathematics has a nonstandard enrichment. 
The enrichment of the naturals is *N, the hypernaturals. 
*N contains N, as well as illimited hypernaturals. 

Not only objects, but (first order) mathematical state- 
ments have ★-transforms. Consider a sequence {xn)n&i 
of real numbers. It has a ^-transform which is a map 
from *N to *R, and which we continue to denote by 
Xn- The statement "lim„_>oo a^n = o", is equivalent to 
" given fc G N, there is such that Vn S N, n > 
a;„ — a| < 1/fc". N, a and are parameters in this state- 
ment. Since *a = a and *mk = mfe, the ^-transform of 
the statement is Vn e *N, n > ruk => \xn — a| < To 
obtain the ^-transform, just "put stars on everything". 
The point of the operation is that the Transfer Principle 
asserts that statement in the standard universe is true 
if and only if its ★-transform is true in the nonstandard 
universe. In the case at hand, an illimited n will satisfy 
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all the statements as k runs through N, so Xn ~ a for 
illimited n. And conversely, if a;„ « a for illimited n, a is 
the limit of the original sequence. 

In a similar way, one finds the nonstandard characteri- 
zation of continuity for a function / : R — > M. Namely, for 
x G M and y E *M. with y ~ x, f{y) « f{x). As just hap- 
pened with * f{y), stars will sometimes be dropped when 
no ambiguity is possible. If x is standard and y ^ x, 
then X is the standard part of j/, and wc write x ^ °y 
or a; = sty. For a standard point x, the collection of all 
points infinitesimally close to it is called its monad. Con- 
tinuity of / at a; is equivalent to "/ maps the monad of x 
into the monad of f{xy\ These concepts of monad and 
standard part generalize to arbitrary topological spaces. 

For another example, consider the least number prin- 
ciple: 

MB E V{N), 3x e B -.y e B ^ X <y. 

This has ★-transform 

VB e *r{N), 3x e B -.y e B => X <y. 

The one free parameter, 7^(N), the power set of the nat- 
urals, has been replace by *V{H). This is not the power 
set of *N. Instead, it is the set of internal subsets of *N . 
A set is standard if it is the ★-transform of a conventional 
set, for example *N. A set is internal if it is an element 
of a standard set. Standard sets are also internal. Not 
all subsets of *N are internal. For example, *N\N is not. 
Indeed, if n is illimited, so is n — 1. Thus, since the set of 
illimited hypernaturals does not obey the least number 
principle, it must not be internal. It is external. Similar 
reasoning leads to the conclusion that the set of infinites- 
imals in *R, or the monad of any point, is an external 
set. 

The internal/external dichotomy is subtle, but cru- 
cially important. If A is an internal subset of *N con- 
taining all illimited n, then there must be some limited 
TO such that n > m ^ n G A. This consequence of the 
least number principle is called underflow (or underspill) . 
Similarly, an internal subset of "^N containing all suffi- 
ciently large limited numbers contains also some illimited 
numbers {overflow or overspill), and an internal subset 
of *R containing all non-negative infinitesimals contains 
all a; < e for some e > in R. 

How can we recognize internal sets, apart from the 
definition? The Internal Definition Principle gives an 
answer. If A is an internal set and (j){x; A,Bi,..., _B„) is 
a formula with internal parameters A, Bi, . . . , _B„. then 
the subset of A consisting of the elements satisfying 0, 
{x G A : (j){x, A, Bi, . . . , Bn)} is also internal. Sup- 
pose an internal sequence (x„)„g*N is given with the 
property that Xn ~ for all limited n. Then, the set 
{to G "^N : \xm\ < l/jTi} is internal. The parameters in it 
are "^N and (xn), both of which are internal. Thus, that 
set is internal by the Internal Definition Principle, and 
therefore contains all to up to some TV in *N\N, implying 
that Xjn ~ for all m < N . ("Robinson's Lemma") 



Another important principle is Saturation. For an in- 
finite cardinal saturated nonstandard model has 
the property that, if {Aj}j^j is a collection of internal 
sets with an index set of cardinality less than k and the 
finite intersection property, then C\j£jAj ^ 0. Usually 
Hi-saturation is adequate for applications. We assume 
that much, at least. 

The sequence space and its nonstandard enrichment 
is a prototype of some of the spaces which will be used 
later. By definition, ti is the Banach space of sequences 
X : N C with finite norm |ja;|| = X^j^i Then, 
*li consists of sequences *N — > *C with norm \\y\\ = 
J2j£*N\yj\ mapping *£i into *M. If x is standard, (in 
£i) then y is in the monad of x (y « x) if \\y — x\\ = 
J2j£*fi lUj^'^jl ~ 0- Here, the ★ has been dropped on the 
nonstandard extensions of both the standard element x 
and the norm. We seek a more explicit condition ensuring 
y « X. 

y is said to be nearstandard if it is in the monad of 
some standard vector in £i, which must therefore be its 
standard part °y. I claim that y is nearstandard pre- 
cisely when 11?/ 1 1 is limited and y "has infinitesimal tail": 
X]i>n \Vi\ ~ fo'' S'll illimited n. And. in that case the 
standard part of y is given by 

{°y)i = °{yi), for limited i. (15) 

The first condition is obvious. The second follows from 
the fact that a standard x € £i has infinitesimal tail: 
given A: G N, there is A''^ such that n > Nk ^ J2i>n 1^*1 ^ 
1/k. By Transfer, n w +oo satisfies all of these state- 
ments as k varies, and therefore X]i>)i ~ 0. (The 
★ on ^*x^ has been dropped again.) So. y ^ x implies 
y also has infinitesimal tail. All that remains to verify, 
then is that pS]) gives the standard part. Define x G ii 
by Xi = °{yi) for limited i. Then, given any natural fc, the 
set Ak ~ {n : IVi ~ < V^} contains all n G N. 

It is also internal, since it has two standard parameters 
(k and x) and one internal but nonstandard parameter 
(y). Thus, there is an illimited N in A^. But then, 

EiG'-N ~Xi\= J2i<N IVi + Ei>Ar IVi -Xi\<l/k 

because N € Ak and both y and x have infinitesimal tail. 
This being true for every fc, \\y — x|| ~ and x is °y, the 
standard part of y. 

The nonstandard characterizations of open, closed and 
compact subsets of a topological space are simple and 
useful. st^^A denotes the union of the monads of all 
points of A. A is open if st^^^ C *A, it is compact if 
*A C st^^^, and it is closed if st^^A contains all near- 
standard points in *A. 



4. LOWER SEMICONTINUITY AND NEARSTANDARD 
MIXED STATES 

This Section establishes lower scmicontinuity of the en- 
ergy forms introduced in Sec. O both for pure states 
and mixed states. Some of the consequences, particularly 
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Corollaries 14.11 and I4.2[ will be important later. The re- 
sults are somewhat technical, but are a nice illustration 
of nonstandard methods. 

We begin with a nonstandard approach to lower semi- 
continuity oi £ on i^, which is a case of the equivalence of 
closedness and lower semicontinuity of positive quadratic 
forms^. Returning to the quadratic forms (|6l7l8p . regard 
them as initially defined on the space of (antisymmetric) 
C^^-valued infinitely differentiable functions of rapid de- 
crease at infinity (the components correspond to spin in- 
dices). That is, (1 + \z\^y"'d'^f/dz'^ is bounded for aU m 
and multi-indices a. Denote this "seed space" by S; other 
choices work, but this one has especially nice properties 
with respect to Fourier transformation. 

It is very important that the form £'~^ is bounded rel- 
ative to £^ with relative bound less than li^^^^, that is, 
for some < a < 1, a sufficiently large b can be found 
such that 

f^(V') <a£^(V') + &llV'f- (16) 
This means that the non-degenerate inner product 

{x\y)£^£ix,y) + {x\y) (17) 

is equivalent to {x\y)K, defined similarly but with £ re- 
placed by £^. That means that there is some c > 
satisfying c~^{-\-)k < {■\')£ < c{-\-)k- One nonstan- 
dard way to construct the completion of S with respect 
to (-j-) or {■\-)k is as a weak nonstandard hull. This 
abstract construction actually applies to completion of 
any linear space S with respect to a non-degenerate in- 
ner product. Let W be the collection of seminorms 
X i-> for cj) € S, initially defined as real-valued 

on S, but immediately extended to hyperreal seminorms 
on *S. Now define the set of VF-infinitesimal elements 
mtw*S of *S to be those satisfying « for every 

(j> G S, and the M^-finite elements finw^S those satisfy- 
ing |(a::|(/')| <C oo for every (j) G S. The completion of 
S is obtained as the quotient S = faiyy* S/mfw* S , and 
the norm is \\x\\ = sup : cj) e S, ||0|| = 1}. It is 

clear that S can be viewed as a subset of S by iden- 
tifying elements with their equivalence classes. To see 
that S is complete, consider a Cauchy sequence (i„) in 
5, lifting to a sequence (x„) in *S. By overspill and 
countable saturation, the sequence can be extended up 
to G "^N \ N still respecting the Cauchy condition 
m > kn =^ \\xm — II < 1/*^ for each n in N. Thus, xn 
is the sought limit of the sequence in S. 

If S is taken to be the completion of S, then this 
constuction makes it crystal clear is that , the comple- 
tion of S with respect to £^ norm can be identified with a 
subspace of S. For, since {x\4>)k = {x\HQQ(j)), where Hqq 
is the kinetic energy Hamiltonian (well-defined as a map 
S — > S*), ('I-) monads are no larger than {•\-)k monads. 
On the other hand, they are also obviously no smaller 
due to the relative weakness of the || • || norm compared 
to the \\-\\k norm. 



On i], the ordinary and £^ norms can be found ac- 
cording to ||a;|| = f{x), \\x\\k = fxix), where f{x) : 
= sup{|((/)|a:)| : ^ G 5,11011 = 1}, and fxix) : = 
sup{|(0|.t)/<-| : (j) £ S, ||(/)||a' = 1}. These take values in 
[0,oo]. Since each seminorm x t-^> \{(j)\x)K\ = \{Hoo4'\x)\ 
is weakly continuous, the supremum is lower scmicon- 
tinuous. Another way to this conclusion is to look at 
the function */x on * fxix) = sn\){\{(j)\x) k\ '■ 4> G 
*S', ||(/)||a- = 1}- Now, if X is standard, that is to say, in 
^, then */(.t) = f{x). For y « x, |(<^|x)k| « |('/>|2/)a| for 
every cj) S, yet the supremum defining \\y\\K includes 
also (j) S *S\S, so that \\y\\K ^ \\x\\k which is exactly the 
nonstandard characterization of lower semicontinuity. 

Lemma 4.1. £i is lower semicontinuous on 9) 

Proof. The proof is contained in the preceding two para- 
graphs. Also, see Reference^. □ 

Turning to mixed states, F G *S is nearstandard if 
it is nearstandard with respect to trace-norm (the only 
topology we consider on ©). It will be useful to charac- 
terize this more concretely in terms of nearstandardness 
in *S^. Whenever possible and profitable we try to reduce 
properties of mixed states to corresponding, but more in- 
tuitively graspable properties of pure states. Now, note 
that for 7 = X^jgn Ci|V'i)('/'i| in given e > in M, there 
is a finite-rank mixed state within e of 7 in trace norm. 
Indeed, some truncation 'Y^^iLiCi\ij}i){'il}i\ will serve the 
purpose. Otherwise put, the finite-rank mixed states are 
dense in ©. Thus, a mixed state 7 = '^^^^:■^c^\'^|)i){'^pi\ 
in *6 is nearstandard if and only if it can be approx- 
imated to any given standard accuracy by some par- 
tial sum Ei6j°Cj|°Vi)(°V'i| with |J| < co. With the 
terms of the sum ordered by decreasing eigenvalue, this 
is equivalent to: ipi is nearstandard whenever Ci 96 0, 
and X^i^i °7 = °Tr7. Also, note that trace-norm topol- 
ogy is stronger than uniform topology, which in turn is 
stronger than strong-operator topology, so that 7 w 7' 
only if every eigenvector of 7 is infinitesimally close to 
being an eigcnector of 7' with the same eigenvalue: if 
(7 — a)V' = 0, then (7' — a)^ ~ 0. 

The way is now prepared for the analog of Lemma 14.11 
for mixed states. 

Lemma 4.2. £1 is lower semicontinuous on & 

Proof. Take 7 G 6, with nonstandard enrichment *7 = 
EiG*NC»lV'i)('/'i|> and suppose 7' = E»e*N CilV'O('0il sat- 
isfies 7' « *7. We need to show that £^1(7') > ^1(7)- 
As just discussed, every eigenvector tp'^ of 7' is infinitely 
close to an eigenvector ipi of 7 with infinitely close eigen- 
value, and vice versa. But, since the eigenvectors of 7 
are standard, Lemma |4. II implies that 

£M) > £iCijl) = £i{ij.). 

(We assume that the terms are arranged in decreasing 
order of eigenvectors and in case of equal eigenvalues. 
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appropriate choices of basis must be made.) Thus, 

m m 
i=l i=l 

for any m £ N. From this, we immediately conclude that 
fi(7')>^i(7)- □ 

The map Dens from to is continuous, therefore if 
-0 G *iO is nearstandard, DensV' is also L^-nearstandard. 
Going the other way, a weaker condition on Densip, to- 
gether with encrgy-boundcdness of ip suffices to guaran- 
tee nearstandardness of ^E* as we shall see. 

If for any e > 0, there is some R such that ^^^^ pdx < 
e, we will say that p is "nearstandardly concentrated" . 
Note that, by underspill this is equivalent to ^^^j^ * p dx ~ 
for all illimitcd R. 

Lemma 4.3. If'^£*Sj has nearstandardly concentrated 
density and is limited, then ^ is nearstandard. 

Proof. Note that, since is a complete metric space, 
it suffices to show that vj/ can be approximated to any 
standard accuracy by a vector in Sj. Also, by relative 
form boundedness of £ with respect to £^ , together 
with positivity of £'-'', there is some constant c such that 
c~^£i < £^ < c£i, so that £ and £^ are interchangeable 
for our purposes here. 

The first step is to show that ^ can be approximated 
by a *Sj vector with bounded support and finite energy 
(a "smooth truncation" ) . Then we will be able to assume 
that itself is such a vector with no loss of generality. 
So, let xi^) be a smooth cutoff function which is nonin- 
creasing with equal to one in Bi and supported in 
B2, where Bji is the ball ||a:|| < R. For R > 0, define the 
iV-particle cutoff function by 

N 

and cutoff wavefunction by 

Thus, has zero probability for any particle to be 
outside the 2i?-ball. By the assumption of nearstan- 
dardly concentrated density, — ^'|| can be made 
smaller than any standard tolerance by taking R large 
enough. Further, since V$fl = XR^"^ + (^X-r)* and 
Xr is smooth, £^{^ii) is limited. Thus, we may assume 
that has bounded support, contained in the closed box 
= {\x\, \y\, \z\ < L}, and wc now consider that case. 
Using the particle-in-a-box wavefunctions for A^, wc 
can construct an explicit orthonormal basis for \l, the 
closed subspace of consisting of vectors with density 
supported in A^. For any E, there are a finite number 
of basis vectors with £^ < E. Thus, given any E and 
e, there is a finite orthonormal set (/ji , . . . , such that 
whenever 4> € \]^ satisfies £i{4') < E, then the compo- 
nent of (j) orthogonal to span((y5i, . . . , ipn) has norm less 



than e. By Transfer, this continues to hold for (/) G *i5 \l- 
But, the projection of (p onto span((^i, . . . , ipn) is clearly 
nearstandard. Since E and e are free, this means that 
can be approximated to any standard tolerance by a near- 
standard vector, and therefore is itself nearstandard. □ 

Here is the central result of this Section. 

Proposition 4.1. //, for T G *S, DcnsF is nearstan- 
dardly concentrated and £i{T) is limited, then F is near- 
standard. 

Proof. Let F have eigenfunction expansion 

r = ^ Ct\ipt){ipt\. 

Since Dens and £1 are additive and positive maps, when- 
ever Ci is not infinitesimal, -ipi is nearstandardly con- 
centrated and of limited energy, hence nearstandard by 
Lemma 14.31 Thus, with 

00 

F':=^°c,|>,)(>,|, 

it follows that °F = F' unless the sum of infinitesimal 
weights Ci is noninfinitesimal, that is TrF' < *TrF — e, 
for some e > in M. In that case, defining J„ : = {i : 
Ci < 1/n}, the internal set K := {n G *N : J^j > e} 
contains all of N, hence some illimited w, and | J^^j > eoj. 
Now define the mixed state 

f ^c,|^,)(V',|<F. 

Since Sj is separable, fewer than ecj/2 of the (orthogonal) 
eigenvectors ipi of F can be nearstandard, and they carry 
a combined weight not exceeding (ew/2)(l/a;) = e/2. 
The remaining eigenvectors are not nearstandard, so 
each must have either illimited £, or non- nearstandardly 
concentrated density, or both. But since these non- 
nearstandard vectors enter F with a combined weight of 
at least e/2, F itself and a fortiori F must also have cither 
illimited £ or non-nearstandardly concentrated density. 
But, that is contrary to hypothesis. □ 

This Proposition has a couple of easy but important 
corollaries. 

Corollary 4.1. If p G X satisfies F[p] < 00, then there 
exists 7 G 6 with Dens 7 = p and f (7) ~ F[p]. 

Proof. Since F[p] < 00, there is some 7 G *6 with 
Dens 7 = p and ^(7) ~ F[p]- By Prop. 14. 1[ 7 is 
nearstandard. Then, since £ is lower semicontinuous, 
bu Lemma [4.21 £{stj) < F[p]. Since both sides of this 
inequality are standard, 

£:(st7) <F[p]. 
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On the other hand, Dens : 6 — >■ X is continuous, which 
imphes that p is in the monad of Dens st 7. But, p is stan- 
dard, so Densst7 = p. Therefore, 8{st"i) > F[p], which 
combined with the previous display yields. £{st"/) ~ 
F[p]. □ 

Corollary 4.2. LetT be a topology on Range Dens which 
makes the map Dens from mixed-states to J^-densities 
continuous, then F is lower semicontinuous with respect 
to T. 

Proof. Let p be a 4-density in Range Dens and p' G 
*RangeDens = Range *Dens such that p' « p. And, 
let 7' be a *density matrix with Dens 7' = p'. By Prop. 
14.11 7' is nearstandard. Thus, Dens "7' is a standard den- 
sity which is infinitely close to p' by the assumption of 
continuity of Dens. Therefore, Dens "7' = ° p' = p, so 
that 

^^[p']=^^i(7')>^i(V)>^[p], 

where the almost- inequality follows from Lemma 14.21 and 
the final inequality by definition of F[p]. □ 

5. COARSE-GRAINING 

This section reviews the basic notions of the coarse- 
grained DFT framework of References^Sr— , augmented 
to the SDFT setting. Much of it looks exactly the same 
on the surface as the conventional fine-grained theory of 
Sec. [21 Building on the results of Sec. HI the existence 
and regularity results for the coarse-grained theory will 
be presented in Sec. |6l In addition to the considerably 
simpler and more intuitive proofs, those results go be- 
yond the previous version of the theory^Sr^^^ by treating 
spin-densities. 

Consider the partition *p of R"^ into a regular array 
of cubical cells fti of side length £. A coarse-grained 4- 
density is simply a specification of the average 4-density 
in each cell. As such, a coarse-grained 4-density is an 
equivalence class of densities in X, the members of the 
class differing only in the way particle number and spin 
is distributed within the cells. Generally, coarse-grained 
densities are denoted by the same sorts of symbols as 
were used previously for fine-grained densities, but there 
should be no confusion as the rest of the paper focusses 
on coarse-grained densities. Since p assigns an aver- 
age density to a cell, it is convenient to identify p with 
the leveled-out function which is equal to that average 
throughout the cell. That makes p into a piecewise con- 
stant function. This is just a very convenient represen- 
tation; the lack of smoothness has no real significance. 
What is done with the (shared) faces of the cubes does 
not matter since they have zero Lebesgue measure. The 
coarse-grained densities now belong to the space of 
cell-constant 4-component functions with norm 

/ (l/o(^)l + |/>)|) dx. (18) 



This Banach space, which will be denoted X, is a sub- 
space oi X, so we could also describe matters with the 
aid of the projection n : X X which averages over 
cells of *p. Then, we are identifying p with 7rp, since irr] 
is the same for every t] ^ p. 

Some useful subspaces of X are singled out: 

X+ -.^ {p e X : po > 0, and \p\ < po} 
X++ := {p e X+ : po > 0} (19) 

X~^ consists of true densities, X^'^ adds the restriction 
that the number-density does not vanish on any cell, 
which is a physically untroubling restriction, as presum- 
ably an infinite potential would be needed to drive the 
density to zero on an entire cell. If jpl = po on some cell, 
that cell is said to be spin- saturated. Subspaces without 
spin saturation 

X® := {p e X+ : \p\ < p„} 
X'^'^ -.^ {p e X++ : \p\ < po} (20) 

are also defined. Only X®® can be fully controlled in a 
general way, but collinear spin-saturated states will be 
treated in Sec. [51 

A fine-grained density in the conventional theory has 
associated with it a whole set of states, and some are 
selected out by an energetic criterion. We do the same 
here. The Lieb internal energy functional is defined as 

F[p] := inf {£:(r) : F G 6, ttDcusF = p} , (21) 

for ph\ X , in perfect analogy to the conventional case. If 
there are no states with ttDcusF = p, then F[p\ — +00, 
as usual. But, there are states realizing any density in 
X'^ , and 

F[p]<\\p\\iVo, forpeA'+, (22) 

where the constant Vb depends on This is in stark 
contrast to the continuum situation, and arises because 
there is a finite intrinsic energy cost A^Vg to putting all 
the particles in the same cell of *p. Any density in X'^ can 
at least be realized by a mixed-state sum of such single- 
cell states, and the bound follows. Since for t G [0,1], 
the convex combination tF + (1 — t)T' of density matrices 
gives the corresponding convex combination of densities, 
it is immediate that F is convex: 

F[tp + (1 - t)p'] < tF[p] + (1 - t)F[p']. (23) 

Now we turn to potentials. For p Cz X, any cell- 
constant 4- vector function v which is bounded below and 
satisfies 

F[p'] + {v,p')>F[p] + {v,p)=0 yp'eX, (24) 

is called a representing potential for p. The reason for the 
name is clear: there is a state F with TrDensF = p with 
total energy in presence of v, and no state with lower 
energy. The stipulation that the ground state energy be 
zero merely fixes the constant which could otherwise be 
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added to vq. It has the useful effect of making the varia- 
tional principle (j24p valid also for improperly normalized 
densities. It is allowed that {v, p') = +00, just as F takes 
value +00 for some arguments. And that may very well 
happen if v is unbounded above, though not for a ground 
state density of v. The set of representing potentials for 
p is denoted by Potl (p). If Potl (p) is non-empty, then p 
is said to be V-representable, abbreviated "VREP" . The 
next Section will show that all of X®® is VREP. 

Because the cells of Vji are non-infinitesimal, there is a 
state of finite intrinsic energy NVq with all particles in 
Hi. Thus, if Do — Iw] < Vq in O^, then that state would 
have total energy less than zero. Thus, our convention 
on the constant offset of potentials implies a floor: 

vo - \v\ > -Vo. (25) 

This can also be written simply v > —Vq if the ordering is 
understood as the usual ordering for Hermitian matrices. 
The set of four-potentials respecting this lower bound, 
which is all that will be needed, is denoted V. Necessarily, 
Potl (p) C V. Similarly, an upper bound can be obtained 
if we also consider the density, li v € Potl (p), then it is 
easy to see that 

K - \v\)po < N^. (26) 

The only topology on V that we will use is the product 
topology. Neighborhoods correspond to constraints on 
the potential on a finite number of cells, and a basis is 
given by sets of the form 

Ue.L{w) ^{veV:\w-v\<e, V|x| < L} (27) 

parametrized by e,L > and to G V, where |k; — t;| = 
wq — vo + Iw — vI. Convergence of a sequence of potentials 
with respect to the product topology is cell-wise conver- 
gence. 

6. V-REPRESENTABILITY AND REGULARITY 

This section is the core of the paper. In it, all the 
basic existence and regularity results are given. In the 
following section, the uniqueness question is taken up. 

A. baby steps 

In Sec. [21 we noted that cpiF = 

{{p\y) & X xR:y > F[p']} is convex and closed, 
but despite that, a straightforward Hahn-Banach ar- 
gument is not available to prove V-representability 
of even one density because domi^ has empty inte- 
rior. This problem is common to the conventional 
and coarse-grained theories. However, in the latter, a 
restriction to a fixed finite volume allows us to make 
the argument work. The conclusion is correspondingly 



weaker, of course, but the simplified problem will serve 
as a foundation on which to build a solution to the 
general problem. 

Consider the densities in which are supported in 
the cube A^, (|a;*| < L for i = 1,2,3) and call this set 
X^ . Then, X^ is a convex set with nonempty interior, 
viewed as a subset of R"'!'^^ I ; this interior consists of den- 
sities which are strictly positive over Aj^ and nowhere 
spin-saturated. The restriction of F to Xj^ is a convex 
function with an upper bound given by inequality (j22p . 
In this case, all the conditions are satisfied to apply the 
separation theorem, and the boundedness of F implies 
that a separating hyperplane at an interior point p is not 
vertical. Thus, it can be represented as a potential v in 
this restricted sense: 

F[p'] + {v, p') > F[p] + {v, p) = Vp' G X+. 

In this case, we say that p is L-VREP. The nomenclature 
is meant to emphasize the property it has which is akin 
to V-representability, even though it is just equivalent 
to being supported essentially on A^ but nowhere spin- 
saturated on that set. 



B. hyperfinite bootstrap 

Using the modest result of the previous subsection, we 
now launch into the hyperfinite. From there, the infinite 
volume limit is an easy step down. Since the L-VREP 
property is standard and holds for every L G N, it holds 
in the nonstandard extension for every L G *N, by the 
Transfer Principle. This simple observation paves the 
way for the following key Lemma. The potentials which 
occur in the hypothesis are in ns*V, that is, in *V and 
near-standard. Since V has the product topology, v G 
nsV means that \v\ > —Vq and v{x) is limited at all 
limited locations x. 

Lemma 6.1. If p £ *X is nearstandard, and is either 
VREP or L-VREP for L ?a +00, with a representing 
potential v G ns *V, then 

(i) . sti; G Potl(stp) 

(ii) . F[stp] « F[stp] 

(iii) . 17 • p is nearstandard in L^, hence v ■ p ?» °v ■° p. 
Proof. According to the hypothesis, 

F[p] + (t;,p) =0, (28) 
for some v G *V. We will show below that 

F[p'] + {°v,p')>Q, Vp'GA'+, (29) 

and that 

{°v,°p)<{v,p). (30) 

With F[°p\ < F[p] from Cor. this yields F[°p] + 
{°V:°p) < 0. Since F[°p] and {°v,°p) are standard, all 
these near-inequalities imply 

F[°p] + {°v,°p)^0, (31) 
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as well as 

F[V]«F[p], {°v,°p)^{v,p). (32) 

This last display gives part (m), and (|29l3ip together 
prove part (i). 

So, to fill in the gaps, suppose that p' violates condition 
(EH): 

F[p'] + [ °v- p'dx<0. (33) 

Then, there is a smooth truncation of p' that also 
violates condition (|29|) . ("Smooth truncation" here refers 
to the construction in Lemma |4.3[ which is really applied 
to the state underlying p', not to the density itself.) For, 
by taking the truncation radius R big enough, F[p'j^] can 
be made as close as desired to At the same time, 

J°v-p'dx = J {°v ■ p')+ dx- J {°v ■ p'Y dx, 

and taking R large makes the change in the integral over 
the negative part as small as needed, due to condition 
(PS)) . while the integral over the positive part can only 
decrease, which is not a threat to the inequality. 

Thus, p' in ([55)) may be assumed to have bounded 
support. But in that case {°v ■ p') sa {v ■ p'), since both 
integrals are effectively only over the i?-ball Br, and v 
is limited there, indeed uniformly bounded both above 
and below on B^- But then F\p'\ + ■ p') < 0. This 
is contrary to the VREP hypothesis on p In case p is 
only L-VREP, the fact that p' has bounded support also 
enters the argument. Thus, condition (|29p is proven. 

Similar reasoning using the lower bound on v and near- 
standard concentration of p proves inequality (|30p . 

(iii): Both p and v are limited at limited points. Thus, 
if p • t) is nearstandard in L^, then °i; • °p is its standard 
part. To show p v is nearstandard, it will suffice to verify 
that 1''^ ' P\dx ~ for illimited L. (The last few 

paragraphs of Sec. [3]are relevant to these considerations.) 

Since p is nearstandard, the lower bound (j25p implies 
that ' P)^ dx « 0. On the other hand, combining 

that with j^^i^^ \*{°v ■ °p)\ « and {°v,°p) w {v,p) 
yields /|^l^^(?; • p)+ rfx « 0. □ 

A general VREP result follows immediately. Recall the 
definition [20] of A"®®. 

Proposition 6.1. Every p e X®® is VREP. 

Proof. Truncate *p, by setting it to zero outside A^ for 
illimited L. This density is L-VREP, and Lemma 16.11 
then gives the result. □ 

Now, by the Transfer Principle, all densities in *X®® 
are VREP, so that Lemma 16.11 applies to any density in 
This observation yields two additional corollaries. 

Proposition 6.2. p i—> F[p] is continuous. 



Proof. Take p a standard density and p' w p. According 
to the transferred version of Prop. 16.11 p' is VREP, and 
by Lemma [ID (m), F[p'] w F[p\. □ 

The second corollary is about continuity-like prop- 
erties of Potl as a function of density. Since p H> 
Potl (p) is a set- valued map, we need to discuss what 
continuity means in that case. A set-valued map / : 
X 1-^ Y between topological spaces is said to be up- 
per semicontinuou^^— aX x & X if, for any open set 
U containing /(x), there is a neighborhood V oi x such 
that f{z) C U for every z ^ V . Upper semicontinuity 
of p H> {v ■ p : V ^ Potl (p)} will be shown to hold on 
X®® . This cannot hold on all of X because Potl takes 
the value at some densities with po = somewhere. 
However, we will show that the graph of Potl is closed 
in A" X V. This means that if (p„, v„) £ GraphPotl is a 
sequence of density/potential pairs with p„ p ^ X m 
and t;„ — ?> £ V cell- wise, then v E Potl(p). 

Proposition 6.3. (a) For every p E X, Potl(p) is 
closed in product topology. 

(b) GraphPotl is closed in X x V. 

(c) The set-valued map p t— > {v ■ p . v E Potl(p)} is 
upper semicontinuous on X®® . 

Proof. Recall the nonstandard characterization of a 
closed set: A is closed if and only if all nearstandard 
points of *A are infinitely close to A. 

(a) Let V G *Potl (p) be nearstandard. Then °v £ 
Potl(p). 

(b) . Let {p,v) be a nearstandard point in 
*GraphPotl = Graph *Potl. This means that both p 
and V are nearstandard. By Lemma l6.1f i). °{p,v) = 
i°p,°v) £ GraphPotl. 

(c) . Let p £ *X®® and v £ Potl(p) both be near- 
standard. Then, p is VREP by the transferred version of 
Prop. 16.11 Thus, by Lemma WAM iii). st (v-p) is infinitely 
close to some element of °p • Potl (°p). □ 

Finally, we show that Potl is the functional derivative 
on X®® , in a sense appropriate to general convex func- 
tionals. Recall that for p, p -\- 5p G X®®, convexity of F 
guarantees existence of the directional derivative 

F'[p;Sp] ■.= \im-{F[p + eSp]-F[p]). 

eio e 

By the very defintion of Potl (p), 

{-v,dp) <F'[p;Sp], 

for every v in Potl (p). What the next result says is that 
this inequality is saturated: 

F'[p; 5p] = max {{-v, Sp) .v e Potl (p)} . (34) 

Thus, all directional derivatives can be recovered from 
Potl (p). If Potl (p) is a singleton, it is easy to see that the 
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potential can be obtained from the directional derivatives 
in a local way: 

dF 

That this remains the case in the non-unique case will 
follow from Prop. 17.21 

Proposition 6.4. For p and p' = p + Sp in X®® , there 
is V Cz Potl (p) satisfying 



Thus, 



v,5p)=F'[p-5p] 



F'[p; 5p] = max {{~v, 5p) -.v £ Potl (p)} . 



Proof. Let e > be an arbitrary infinitesimal. Then, 
since (one-sided) directional derivatives of F exist, 

F[p + e8p\-F[p\=e{F'[p-5p]+T^), 

where ~ 0. Now, p + eSp is VREP by Prop. 16.11 and 
Transfer, so let v S Potl (p -I- eSp). Then, 

F[p + edp]-{v,-edp)<F[p]. 

Together with the previous displayed equation, this yields 

{-v,dp)>F'[p;Sp]. (35) 

On the other hand, by Lemma [6.11 °v G Potl(p), im- 
plying that 



{-°v,5p) < F'[p;Sp]. 



(36) 



We would like to get equality here by showing that the 
near-inequality psp is, roughly speaking, nearstandardly 
carried, in order to replace "<" in with ">" . 
Suppose then that this hope is dashed and 



V ■ 6pdx <^ 0, 



(37) 



x\>L 



for some L ^ oo. By the argument in Lemma |6. II 



V ■ {p + eSp) dx K, 0, 



x\>L 



otherwise truncation would show that p + e6p is not the 
ground state density of v. But, then, by combining the 
last two displays, 



V ■ {p + Sp) dx < 0. 



(38) 



x\>L 



Now, a "backward truncation" of p + Sp, keeping just 
the part for \x\ > L would have infinitesimal intrinsic en- 
ergy, but potential energy psp appreciably less than zero, 
again implying that p-\-e5p is not a ground state density 
in V. Thus, (|57)) is false, implying that the inequality 
psp is nearstandardly carried and the inequality in (|36p 
can be turned around. This proves the Proposition. □ 



7. FOUR-POTENTIAL UNIQUENESS AND 
NON-UNIQUENESS 

This section is concerned with the question of when 
Potl(p) might have more than one element. This 
problem attracted renewed scrutiny in the conventional 
settingi^"— , long after the original investigation^. The 
conclusion reached here (Prop. 17. 2p is similar to the cu- 
mulative conclusion of those works. But here unique 
continuation properties are demonstrated rather than 
being assumed. In quick summary, disallowing spin- 
saturation leaves the only possibility of non-uniqueness 
being collinear magnetization in an Sz eigenstate. The 
next section looks more closely at the collinear case with 
spin-saturation. No nonstandard analysis is used in ei- 
ther of these two sections. 

In contrast to previous sections, it is necessary to care- 
fully consider the behavior of wavefunctions in configura- 
tion space. Throughout this section, is a fixed ground 
state wavefunction of total energy in the presence of 
potential v = (it, B): 



{Ho + y)* = 0, 



where 



V 



N 

E 



^u{Xa) + B{Xa) 



(39) 



(40) 



Here, contains kinetic energy and interaction energy 
and is diagonal in spins. The potential at the posi- 
tion of particle a is {u{xa),B{xa))- 4* might be a zero- 
energy ground state for some other potential. Elements 
of Potl {^) other than (u, B) are generically denoted as 
{u',B'). It will soon become clear that, in studying the 
non-uniqueness problem, only the magnetic field com- 
ponent really requires attention. To facilitate that, the 
notation 

Potl (p) = : (u, B) G Potl (p) for some w| (41) 

will be used for the projection of Potl(p). Naturally 
extending previous notation, Potl (4*), denotes the set of 
four-potentials having ^I* as a zero-energy ground state. 
As previously, position space cells ( "*P-cells" ) are denoted 
fli, etc., i being simply an abstract cell index. Cells in 
configuration space ( "Cp^-cells" ) are denoted in the style 
rif^...j^. In cell particle 1 is in . . ., particle 

N is in ilij^ . The value of {u, B) in cell fi^ is denoted 
{ui,Bi). In this section and the next, instead of (po,p) 
for the 4-density TrDens 5*, we use the notation {ni, 2mi) 
(Note the factor of 2). If for every {u',B'^) G Pot^^-), 
B'l = Bi, then z is a unique- cell (UB cell), otherwise a 
non-uniquc-S cell (non-UB cell). 

The following unique continuation result for single- 
component wavefunctions is fundamental. 
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Thm. XIII.63 of RSIVi^. Let u G Hf^^, that is, 
(j)u e D{-A) for each (f> £ C^(R"). Let D he an open 
connected set in M" and suppose that 



\V'^u{x)\ < M\u{x)\ 



(42) 



Now take the commutator of Hq + V with the new oper- 
ator, to find 



TV 



almost everywhere in D. Then, ifu vanishes in the neigh- 
borhood of a single point xq € D, u is identically zero in 
D. 

Al is supposed to be an arbitrary but fixed constant. 
An examination of the proof of the basic unique continu- 
ation theorem shows that it extends to multicomponent 
wavefunctions, such as wavefunctions with spin indices, 
by replacing [up and |V^up by the corresponding sums 
of squares of spin components. Indeed, this replacement 
can simply be made globally throughout the proof of the 
theorem and lemmas leading up to it. Unfortunately, 
that does not immediately give us any control of indi- 
vidual spin components. The usual way a bound of the 
sort (|42|) arises, of course, is that u satisfies a Schrodinger 
equation, so that M is \V — E\. An individual spin com- 
ponent might satisfy such an equation by itself. If it is 
appropriately related to a different component that does 
so, that can also work. 

It may appear at first that the unboundedness of the 
Coulomb repulsion could cause problems. On any A^, 
the single-particle external potential is bounded, so the 
total external potential is bounded on A^. The sum of 
the external potential and the interaction is then uni- 
formly bounded on = A^ n [/„ for any n, where 
Un = {x : \xi - x.j\ > l/n,Wi ^ j}. Also, UL,n is a 
connected (open) set for large enough n, so the theorem 
will apply to it. But, any open subset of A^ must have 
an open intersection with for some n. So, in fact, 
we are free to take D to be A^. This depends on the 
set where the interaction energy is infinite having small 
codimension and no interior. Our first application of the 
theorem is then to note that if Vlf = on any r2f^...j^ , then 
it vanishes everywhere. Thus, there is nonzero probabil- 
ity to be in any *P^-cell. 

Note in passing that, by basic elliptic regularity 
theory^, 'I' is continuous away from coincident configu- 
ration points {xa = xp for a ^ (5) and is C°° away from 
coincidence or particles on the boundaries of *!p-cells. 

If {u',B') € Pot^^-) is distinct from {u, B) then with 
Au = u' - It and AB ^ B' - B, 



N 

Af * = [Au(.t„) + AB{xa) 



= 0. 



(43) 



It then follows that [Hq + V, AV]^ = 0. As a distribu- 
tional equation on the interior of the configuration space 



cell il, 



this yields 



JV 

E 

a=l 



{B,^ X AS, J ■ 0. 



These equations are useful in the case of VLn...i, Our basic 
unique continuation result assures that does not vanish 
on any open set in this cell, so the equations are not 
vacuous. 



First, note from Eq. (j43l) that if AB^ = 0, then Au; = 
as well. So a UB cell also has unique u, explaining the 
focus on the UB/non-UB dichotomy. Thus, suppose that 
Vti is a non-UB cell. Then, ABi can be taken nonzero, 
and so can Bi (otherwise fiip the roles of Bi and B'^). If 
ABi is not parallel to Bi, that makes ^' an eigenstate of 
spin along two perpendicular axes, which is impossible. 
Thus, ABi is parallel to Bi. 

Lemma 7.1. // is a non-UB cell, then there is a 
unique axis ii, such that for every B in Pot^^*), Bi is 
parallel to Ci, that is, Bi = BiCi. 



Proof. See previous paragraph. 



□ 



If is a UB cell with Bj ^ 0, there is also a unique 
axis Cj, namely B/\B\. If = is unique, Cj is arbi- 
trary. Each cell flj thus has its own spin quantization axis 
Cj. In the following, we will mostly use this ei-&aszs. That 

(si • • • SArj^I') has 
= +1 and down 



is, in the component = 

the spin of particle a up along e^^ if s^ 
along ii^ if s^ — —1. 

Consider now the restriction of ^E* to ^il...i,^- The ex- 
ternal four-potential is diagonalized in the e basis, so the 
unique continuation theorem quoted above can be ap- 
plied separately to each component ^'s^...s„ to see that if 
^'sj...s„ vanishes on any open set in rif^...^^, it vanishes 



N 



almost everywhere on the cell Q -^ 

At this point, we want to consider fi^.. 
ik+i, . . ■ ,iN all different from j 



nu^^■■■^N with 



Lemma 7.2. 



If rij is non-UB, then, for fixed 
with ik+i, ■ ■ ■ ,iN 0,11 different from j, the 
Sk+i . . . sn component of"^ (in e basis) is an eigenstate of 
spin along ij. That is, for only one si, . . . , Sk is ^si - sjv 
nonzero. 

Proof. This follows straightforwardly from the eigenfunc- 
tion Eq. ABj can be assumed to be nonzero along 

ej, giving 

N 

kAuj + ABj{k - 2Nji) + ^ (Aui^ + AB,^Sa) = 0, 

a=k+l 

where Nji is the number of spin-down particles in flj. 
The equation has at most one solution for Nj^. □ 

Proposition 7.1. Ifilj is a non-UB cell, then Cj xrhj ~ 
0. 



13 



Proof. This follows from Lemma 17.21 by summation over 
cell occupations and integration. □ 

At this point we turn attention from a single config- 
uration space cell to a pair of neighboring cells ^^^^....^j^ 
and ^ji^-.-ij^ ■ "Neighboring" means that flk shares a face 
with Qj. The indices 12 • • • *Ar are fairly inert, so the cell 
index notation is abbreviated to i^^jj and f2^^j . 

For the following argument we use a slight modification 
of the e basis, namely using in both Qk and flj. {ek[i]- 
basis) . 

Lemma 7.3. Suppose that flj and Qk o,re neighboring 
cells. Using ekii\-basis in both fJ^^j and rij^Jjj, if'i'si - SN 
vanishes in fl^, ], Then it also vanishes on fl^r-\- 

Proof. The unique continuation theorem is applicable be- 
cause the potential is diagonalized in fi^^j , while the con- 
dition p2|) for our chosen component is satisfied in 
by hypothesis. □ 

A possibly useful mnemonic is to say that "moving one 
particle from one cell to another cannot flip the spins of 
other particles" , bearing in mind that "moving" does not 
refer to a dynamical process but just a shift of attention. 

Lemma 7.4. // flj and flk are neighboring ^-cells and 
both non-UB, then Sj and are parallel. 

Proof. In the e spin basis, find an occupied s[t] spin state 
in fijl^j , with all of ^2, . . . , different from j and k. Then 
(— s)[i] is unoccupied by Lemma [7T2] Also, there is some 
occupied s'[t] in f^^^j, in e^-^i] basis, for if there were 
no such s', then Lemma 17.31 would guarantee that s[t] is 
unoccupied in f^^jj- Furthermore, (— s')[f] is then unoc- 
cupied in il^^j , again by Lemma 17.21 
Thus, in n^^^, 

= (-sls')*,,!^] + (-s| - s')*(-.')[t] 
= (-s|s')*.'[t]. 

^'(_s)[t] = in f^^j], and therefore satisfies a bound of 
type (|42|) there. But, it also satisfies such a bound in 
ri^jj. This follows because ^'s'[t] satisfies such a bound 
there by virtue of obeying a Schrodinger equation, and 
^'(_s)[t] has a fixed proportionality to 5's'[t] in fij^jj ac- 
cording to the previous display. 

The basic unique continuation theorem then says that 
^'(_s)[t] = in J^^j]. But that implies that \s) ~ \s'), and 
since these spin states are eigenstates of ej ■ a and ik ■ <J 
respectively, Cj and are coaxial. □ 

The picture according to what has been proven so far 
is of connected components of non-UB cells, separated by 
UB cells. On each connected cluster, there is a unique 
common e axis. So, it is now time to consider UB cells. 
It turns out that the existence of even one has drastic 
consequences. 



Lemma 7.5. If there is a UB cell, then all non-UB cells 
are spin- saturated. 

Proof. Label one of the unique- i? cells '0'. Then, con- 
sidering the *p^-cell i^^..o with all particles in JIq, the 
eigenfunction Eq. (|43| becomes iVAuo = 0. So, Auq = 
also. 

Let rij be a non-UB cell and look in f2j^...g. (If there 
are no non-UB cells, the Lemma is vacuously true.) The 
argument now is like that in Lcmma l7.2l 

Auj + ABj • CTi = 

for any spin component which occurs in the wavefunction 
in fij^...g, where AB^ can be taken non-zero along ej. 
This equation can be satisfied for only one value of s^ = 
ij •(?, either -1-1 or —1 but not both. Furthermore, Auj = 

~s^\ABj\. Shifting attention to ^fjo...o, in 
turn, the same argument shows that all the spins in fij 
arc in the spin state. 

But, this argument docs not tell us anything about 
the wavefunction in Ojfco - o, where fij and arc two 
different non-UB cells. If \ABj\ = \ABk\, it looks as 
though it might be possible to flip a spin in each cell 
while continuing to respect Eq. (|43|) . 

Lemma 17.31 comes to the rescue here. For exam- 
ple, is saturated, by the preceding argument, and 
^f...jk can be reached by moving (shift of attention) par- 
ticle N from one cell to another until k is reached. Ac- 
cording to Lcmma [7.3l at no step in this process do we see 
flips of the particles left in O ^ . so that the only possibility 
allowed by the eigenfunction Eq. (|43)) is that particle N 
ends up in state . It is clear how the example general- 
izes. 

Thus, if there is even one unique-B cell, then any 
particle which is ever in the non-UB cell Q,j is in spin 
state s^ . That shows spin saturation in cell j, namely, 
|mj|=nj/2. □ 

From the foregoing results, the main conclusion of this 
section can now be assembled. If there are UB (unique- 
B) cells as well as non-UB cells, then the latter are spin- 
saturated. We reconsider spin saturated cases in the next 
section, but for now assume no spin saturation, so this 
case is ruled out. All cells are thus non-UB, and all Cj 
are therefore along the same axis by Lemma 17.41 In that 
case, there is a global spin quantization axis which diag- 
onalizes the Hamiltonian. 5^, the total spin along the 
common axis, is a good quantum number. Taking to 
be an cigenstatc of Sz, and looking at fi^.-.i, we see that 
Aui = —{2Sz/N)ABi Thus, different Sz eigenstates can- 
not share more than one potential. Also, using this, the 
eigenfunction Eq. becomes 

which impUcs that all ABi arc equal. Therefore, 
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Proposition 7.2. if p is nowhere spin- saturated but 
Potl (p) is not a singleton, then rn is everywhere along a 
common axis, rn = me, and the total spin J mdx = Sz is 
a half-integer. Furthermore Potl (p) is a one-parameter 
family with AB only allowed to be uniform: 

u, = u1- 2SzAB/N, B, = e[B^ + AB]. 

Proof. Preceding discussion. □ 

Note that this last result implies that Potl (p) can be 
locally determined from directional derivatives for varia- 
tion of the density in only a single cell since such varia- 
tions can be used to determine the ranges of Ui and Bi 
and the proposition shows how to put them together to 
construct all of Potl (p) . 

8. COLLINEAR STATES AND SPIN-SATURATION 



\Sz = ~N/2) is an eigenstate of any of these v. Since Sz 
is a good quantum number, all that needs to be done is to 
add a large enough uniform constant to B (and subtract 
it from u) so that the ground states in all other Sz sectors 
have greater energy than does 7 |S'z — —N/2). 

To see that other cases are not V-representable, as- 
sume m = mz with m = —n/2 somewhere, but not 
everywhere, so that {u, B) S Potl (p) satisfies B = Bz 
according to Prop. 18.11 Sz is a good quantum number, 
so consider the ground state energies in each Sz sector. If 
Sz = —N/2 is not tied for minimum energy, the assumed 
situation is impossible. If a degeneracy does occur, then 
the Sz = —N/2 components can be removed from the 
state and the result is still a ground state for (m, B). But 
this state has zero density in at least one cell, and that 
is impossible because it violates our basic unique contin- 
uation principle. □ 



This short section takes up densities which are some- 
where spin-saturated, especially the important case of 
collinear spin density. As for the general case, it is easy 
to see that a spin-saturated cell is non-UB, so according 
to the remark following Lemma (7.41 the saturated cells 
form clusters all magnetized in the same direction sepa- 
rated by non-spin-saturated, non-UB cells. Prop [6731 fb) 
offers a possible route to find representing potentials for 
such densities. 

The rest of this section concentrates on the collinear 
case. Thus rfi is everywhere along the z axis, though 
it might be zero. Spin saturation has not been allowed 
in the general results of the previous two Sections, but 
collinear spin-saturated states clearly exist in nature, so 
it is important to make special provision for them. 

Proposition 8.1. If rn is everywhere along the z axis, 
then so is B for any B G Potl (p). 

Proof. Assume tt Dens^* = (n, 2m) with rfi everywhere 
along z, and {u,B) G Potl (n, 2m). Simultaneous rota- 
tion of the spins and of B by the same angle about the z 
axis preserves those relationships, as well as £{'i'). Since 
m is invariant under the rotation, if Bt is not, then Qt is 
a non-UB cell. But in that case, B is along z according 
to Prop. O □ 

Proposition 8.2. If p ~ {n, 2m) G and rn is ev- 

erywhere along the z axis, then Potl (p) ^ if and only 
if p is either everywhere spin-saturated or nowhere spin- 
saturated. 

Proof. The case of nowhere spin-saturated p has already 
been dealt with, so suppose that p is everywhere spin- 
saturated with m = —in/2. 

In the spinless version of the theory of Sec. ^ n = po is 
V-rcpresentable; it is the density of some ground state 7 
in potential v. Now, consider the family of four-potentials 
with B = Bz, and u — B = v. The product of 7 and 



9. CONCLUDING REMARKS 

This paper extends previously established^"— results 
of coarse-grained DFT to SDFT. Here is a quick sum- 
mary. As a function of coarse-grained four-density with 
norm, the Lieb internal energy functional F is contin- 
uous. For a coarse-grained density p, Potl (p) is the set of 
representing four-potentials which are constant on cells 
with additive constant fixed by F\p\ -\- {v, p) = 0. Each 
four-potential v in Potl(p) has a ground state 7 with 
a four-density that coarse-grains to p: tt Dens 7 = p. 
If p is everywhere nonzero and nowhere spin-saturated 
(p G X®®), then p is VREP: Potl(p) ^ 0. In fact, 
Potl is the functional derivative of F in this sense: for p 
and p + Jp in X®® , the directional derivative F'[p; 6p\ 
is equal to max {{—v,5p) : v G Potl(p)}. Furthermore, 
p ^ {p ■ V : V € Potl (p)} is an upper semicontinuous 
set- valued function, Potl takes values in closed sets and 
Graph Potl is closed in A" x V. Apart from spin-saturated 
densities, only for densities p with collinear magnetiza- 
tion that integrates to a half-integer can Potl (p) have 
more than one element. And in that case, there is a 
global axis e such that B is everywhere along e for all 
B G Potl (p) and different elements of Potl (p) differ by 
a uniform shift along e. Collinear states with fully satu- 
rated spin are also V-representable. 

This paper has an additional aim of promoting the use 
of nonstandard analysis in mathematical and theoreti- 
cal physics. The results of Section [5] were obtained with 
the aid of nonstandard analysis tools, resulting in proofs 
which are intuitive and of modest technical sophistica- 
tion. The turn to infinitesimal methods was spurred by 
the hope that it would shed light on the limit of coarse- 
graining scale going to zero. In the meantime, it has 
shown its value for the infinite- volume limit. 
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